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ABSTRACT. We outline the physical picture of the high-redshift Lyct forest 
that has emerged from cosmological simulations, discuss statistical character- 
istics of the forest that can be used to test theories of structure formation, 
' present a preliminary comparison between simulation predictions and mea- 

^-j- ' surements from a sample of 28 Keck HIRES QSO spectra, and summarize the 

results of a recent determination of the linear mass power spectrum P(k) at 
2 = 2.5 from a sample of 19 moderate resolution QSO spectra. The physical 
picture is simple if each QSO spectrum is viewed as a continuous non-linear 
qq ' map of the line-of-sight density field rather than a collection of discrete ab- 

, sorption lines. To a good approximation, the relation between Lya optical 

depth and mass overdensity is r = A(p/p) 1,6 . The constant A depends on 
poorly known physical parameters, but for a specified cosmological model its 
O-H value can be calibrated using the mean flux decrement (D) = (1 — e~ T ), 

leaving all other statistical properties of the forest as independent model pre- 
dictions. The distribution of flux decrements is closely tied to the probability 
distribution function (PDF) of the underlying mass distribution, and hence 
to the amplitude and PDF (Gaussian or non-Gaussian) of the primordial 
density fluctuations. The threshold crossing frequency, analogous to the 3-d 
genus curve, responds to the shape and amplitude of the primordial P(k) 
and to the values of f2o and £7a- The predictions of open CDM and ACDM 
models agree well with the measured flux decrement distribution at smooth- 
ing lengths of 25 km s" 1 and 100 km s — 1 and with the threshold crossing 
frequency at 100 km s~ x . Discrepancy with the observed threshold crossing 
frequency at 25 km s — 1 smoothing may reflect the combined effects of noise 
in the data and limited mass resolution of the simulations. The moderate 
resolution spectra yield the slope and amplitude of the linear mass P(k) at 
27r/fc ~ 700 km s" 1 . The slope, never previously measured on these scales, 
agrees with the predictions of inflation+CDM models. Combining the ampli- 
tude with COBE normalization imposes a constraint on these models of the 
form Qoh a nPQ^ ) = constant. Assuming Gaussian primordial fluctuations and 
a power spectrum shape parameter T f» 0.2, consistency of the measured P(k) 
with the observed cluster mass function at z = requires C = 0.461q 12 for 
n A = and C = 0.34+""^| for U A = 1 - Ho (Ict errors). 
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1 Introduction 

For decades, the study of large 
scale structure was virtually synony- 
mous with the study of galaxy clus- 



tering. The 1980s and 1990s have 
seen spectacular improvements in 
galaxy redshift and peculiar veloc- 
ity data, many of them documented 
at this meeting. At this point, our 
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ability to draw cosmological conclu- 
sions from these data is often limited 
less by their statistical uncertainties 
than by our limited understanding of 
the complex physics of galaxy for- 
mation. CMB anisotropy measure- 
ments have been enormously influen- 
tial because their underlying physics 
is simpler, allowing a straightforward 
connection between theory and ob- 
servation. However, even the ambi- 
tious CMB experiments now under- 
way will yield only a particular pro- 
jection of the evolution of cosmic 
structure, mostly concentrated at 
z ~ 1000. Other measures that probe 
3-dimensional structure at later red- 
shifts are needed to complement 
these experiments. 

The success of hydrodynamic cos- 
mological simulations in explaining 
the basic properties of the Lye* for- 
est (Cen et al. 1994; Zhang, Anni- 
nos, &: Norman 1995; Hernquist et 
al. 1996; Wadsley & Bond 1996; The- 
uns et al. 1998) has opened up a 
new testing ground for cosmological 
theories. Ground-based observations 
yield superb data on the Lya forest 
at z ~ 2 — 4, and HST spectra pro- 
vide lower resolution measurements 
at z ~ — 2. Equally important, the 
simulations lead to a simple physical 
picture of the Lyct forest, one that is 
well approximated by semi-analytic 
models that were to a significant ex- 
tent developed prior to and indepen- 
dent of the simulations (McGill 1990; 
Bi 1993; Bi, Ge, & Fang 1995; Bi 
& Davidsen 1997; Hui, Gnedin, & 
Zhang 1997). We will briefly describe 
this physical picture in the next 
section (more extended discussions 
along these lines appear in Wein- 
berg et al. 1997 and Weinberg, Katz, 
& Hernquist 1998, hereafter WKH). 
We will then discuss the physical sig- 
nificance of two statistics that can 
be used to test cosmological models 
against Lyct forest data and present a 
preliminary comparison between pre- 
dictions from smoothed particle hy- 
drodynamics (SPH) simulations and 
measurements from a sample of Keck 
HIRES spectra. We conclude with 
a discussion of a recent determina- 
tion of the mass power spectrum at 
z = 2.5 and its implications for the 
value of Ho and the parameters of 
CDM models. 



2 Physics of the Lya Forest 

In hydrodynamic simulations of 
high redshift structure (z ~ 2 — 5), 
the Lya forest is produced mainly 
by gas with overdensity p/p < 10. 
Most of this gas is unshocked, and 
the competition between photoion- 
ization heating and adiabatic cool- 
ing leads to a tight relation between 
p and T that is well approximated 
by a power law, T = T (p/p) a . The 
values of To and a depend on the 
reionization history of the universe 
and on the spectral shape of the UV 
background; they typically lie in the 
ranges 4000 K < T < 10, 000 K and 
0.3 < a < 0.6 (Hui ^Gnedin 1997). 
The optical depth for Lya absorption 
is proportional to the neutral hydro- 
gen density, which for gas in pho- 
toionization equilibrium near 10 4 K is 
proportional to p 2 T~ - 7 /Tm, where 
Thi is the HI photoionization rate 
and the T -0 - 7 factor accounts for the 
temperature dependence of the re- 
combination rate. The combination 
of photoionization equilibrium and 
the p — T relation therefore leads to 
a power law relation between optical 
depth, t, and gas overdensity p/p, 



T = A(p/pf , (1) 




where (i = 2 - 0.7a RJ 1.6. 

Since equation (|j|) describes the 
analog of Gunn-Peterson (1965) ab- 
sorption for a non-uniform, pho- 
toionized medium, we refer to it as 
the Fluctuating Gunn-Peterson Ap- 
proximation, or FGPA (see Rauch et 
al. 1997; Croft et al. 1998a; WKH). 
The FGPA assumes that all gas lies 
on the temperature-density relation, 
and in the form expressed here it 
ignores the effects of peculiar ve- 
locities and thermal broadening. In 
principle, the quantity p/p in equa- 
tion (hi) represents the gas overden- 
sity. However, in the cool, low den- 
sity regions that are responsible for 
most Lya forest absorption, pres- 
sure gradients are usually small com- 
pared to gravitational forces. The gas 



Evolution of Large Scale Structure / Garching August 1998 



Cosmology with the Lya Forest 



3 



redshift space 



z = 3 



real space 



z = 3 




Figure 1. A test of the Fluctuating Gunn- Peterson Approximation on artifi- 
cial spectra extracted from an SPH simulation of standard CDM in real space 
(right) and redshift space (left), at z = 3. Points show the relation between 
Lya optical depth and gas overdensity in the simulated spectra, pixel by pixel. 
Diagonal lines show the prediction of equation (In). 



and dark matter therefore trace each 
other quite well, with the gas dis- 
tribution being slightly smoother in 
the neighborhood of low overdensity 
peaks of the dark matter distribu- 
tion (Gnedin & Hui 1998; Bryan et 
al. 1998). 

Figure [j] tests the FGPA against 
results from an SPH simulation of 
"standard" CDM (hereafter SCDM, 
with Q = 1, h = 0.5, as = 0.7), at 
2 = 3. The right hand panel plots the 
Lya optical depth against gas over- 
density for spectra extracted in real 
space, i.e., with thermal broadening 
and peculiar velocities set to zero. 
Most points lie on a tight, clearly de- 
fined sequence. Points lying well be- 
low this sequence come from regions 
where shock heating has raised the 
gas temperature above the power law 
p — T relation, depressing the recom- 
bination rate and hence the neutral 
hydrogen optical depth. Peculiar ve- 
locities and thermal broadening in- 
crease the scatter of the relation in 
redshift space (left hand panel), but 
they do not destroy it. The diago- 
nal lines show the prediction of equa- 
tion (Q), with all "free" parameters 
(fi b , h, T , r m , P, H(z)/H ) set 
to the values that they have in the 
simulation. The FGPA provides a 
good but not perfect description of 
the numerical results, breaking down 
mainly in the regions where shock 
heating is important. 

For an approximate numerical ap- 
proach that is much cheaper than a 



high resolution hydrodynamic calcu- 
lation, one can run a lower resolution 
particle-mesh (PM) N-body simula- 
tion, compute the density field from 
the evolved particle distribution, im- 
pose the p — T relation, and extract 
spectra. This technique uses a fully 
non-linear solution for the density 
and velocity fields, but it still as- 
sumes that gas traces dark matter 
and that all gas lies on the p — T re- 
lation, approximations that are good 
but not perfect. As shown in WKH 
(Figure 8), the agreement between 
the PM approximation and a full hy- 
drodynamic simulation is very good 
over most of the spectrum, but it 
breaks down in some higher density, 
shocked regions. Similar techniques 
have been used by Gnedin & Hui 
(1998), who also incorporate an ap- 
proximate treatment of gas pressure 
in the N-body calculation, and by 
Petitjean, Miicket, & Kates (1995) 
and Miicket et al. (1996), who use 
a heuristic method to incorporate 
shock heating. 

3 Continuous Field Statistics 

The approximate physical picture 
just outlined implies a simple and di- 
rect relation between observed flux 
and mass overdensity, 

L = e -r = e-^P/p) 1 - 8 . (2 ) 
F c 
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The unabsorbod continuum flux, F c , 
must be estimated by fitting a con- 
tinuum to regions of low apparent 
absorption or by extrapolating the 
continuum from the red side of the 
Lya emission line. The constant A 
depends on poorly known parame- 
ters, but given a cosmological model 
that predicts the statistical proper- 
ties of the density and peculiar ve- 
locity fields, the value of A can be 
fixed by matching one observable, 
the most obvious choice being the 
mean flux decrement (D) = (1 — 
e~ T ). This mean decrement normal- 
ization calibrates the one unknown 
parameter that determines the rela- 
tion between mass fluctuations and 
optical depth fluctuations. 

The traditional approach to anal- 
ysis of the Lya forest is to decompose 
each spectrum into lines and com- 
pute the statistical distributions of 
the line properties. However, while 
the relation between flux and mass 
density is simple, the relation be- 
tween the number of lines and mass 
density is not. Line decomposition 
also interposes a complicated, non- 
linear algorithm in the path between 
data and model predictions, an al- 
gorithm whose performance depends 
in an intricate way on the signal-to- 
noise ratio and spectral resolution of 
the data. Statistical measures that 
treat the spectrum as a continuous 
field are likely to be more powerful 
than line decomposition methods for 
discriminating cosmological models, 
because they involve minimal manip- 
ulation of the observational data and 
because they are better attuned to 
the underlying physics of the absorb- 
ing medium. 

One of the simplest of these con- 
tinuous field statistics is the distri- 
bution function of flux decrements 
P(D), where D = 1 — e~ T and 
P(D)dD is the fraction of pixels with 
flux decrement in the range D — > 
D + dD (Miralda-Escude et al. 1996, 
1997; Croft et al. 1997; Kim et al. 
1997; Rauch et al. 1997). This statis- 
tic is analogous to the galaxy counts- 
in-cells distribution frequently used 
in studies of large scale structure. 
P(D) can be measured directly from 
unsmoothed spectra, but one can 
also investigate larger scale cluster- 
ing by smoothing the spectra and 
then measuring P(D), just as one 
would measure galaxy count distri- 



butions in cells of different sizes. 
Equation (^) implies that P(D) 
should be closely related to the prob- 
ability distribution function (PDF) 
of the underlying mass fluctuations. 
For models in which the primordial 
fluctuations have a Gaussian PDF, 
as predicted by most versions of in- 
flation, the PDF of the non-linear 
mass fluctuations depends mainly on 
the amplitude of the linear power 
spectrum. Non-Gaussianity of the 
primordial fluctuations should also 
affect the mass PDF, and hence the 
flux decrement distribution. 

Figure || confirms these expec- 
tations. In the upper panels, the 
heavy and light solid lines com- 
pare the flux decrement distribu- 
tions measured from a full SPH sim- 
ulation of the SCDM model and a 
PM simulation with the same ini- 
tial conditions. This comparison in- 
dicates that the PM approximation 
performs quite well for this statis- 
tic at these smoothing lengths (25 
and 100 km s — 1 ). The other lines in 
the upper panels show results from 
PM simulations in which the initial 
mass fluctuations were multiplied by 
a factor of 0.25 (dotted), 0.5 (short- 
dashed), and 2.0 (long-dashed) prior 
to evolution. In all cases the spec- 
tra are normalized to the same 
mean flux decrement. As expected, 
P{D) becomes steadily broader as 
the mass fluctuation amplitude in- 
creases. The lower panels show the 
effect of changing the PDF of the 
primordial fluctuations while keep- 
ing the power spectrum fixed. We 
generate non-Gaussian initial condi- 
tions for the simulations using the 
non-linear transformation method of 
Weinberg & Cole (1992). The skew- 
negative and broad models, both of 
which have extended non-Gaussian 
tails of negative fluctuations, pro- 
duce Lya forest spectra with large 
numbers of low absorption pixels. 
The skew-positive model, by con- 
trast, predicts very few transparent 
regions and a flux decrement distri- 
bution that is more narrowly peaked 
about its mean value. 

Using Eulerian hydrodynamic 
simulations and SPH simulations, re- 
spectively, Bryan et al. (1998) and 
Theuns et al. (1998) have examined 
the effects of mass resolution and 
simulation box size on predictions 
of the unsmoothed flux decrement 
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Figure 2. The influence of mass fluctuation amplitude and primordial PDF 
shape on the flux decrement distribution function, P(D), with smoothing 
lengths of 25 km s" 1 (left) and 100 km s^ 1 (right). In all panels, light solid 
lines show results from a PM simulation of the SCDM model with Gaus- 
sian initial conditions. In the upper panels, heavy solid lines show the results 
from a full SPH simulation with the same initial conditions, and other line 
types show results from PM simulations in which the amplitude of the ini- 
tial conditions has been scaled by the factor a indicated in the legend. Lower 
panels compare SCDM to models with the same initial power spectrum but 
non-Gaussian PDFs that are skew-positive, skew-negative, or symmetric with 
broad tails. 



distribution. Both groups find that 
this quantity is fairly but not com- 
pletely robust to changes in numer- 
ical parameters. Figure |i| presents 
a similar investigation at smoothing 
lengths of 25 and 100 km s —1 , using 
the PM approximation. The speed 
and simplicity of the PM approxima- 
tion make it a useful exploratory tool 
for such a study, although its conver- 
gence properties will not be identical 
to those of full hydrodynamic simu- 
lations because it ignores the influ- 
ence of gas pressure. Figure |5j sug- 
gests that the parameters of our cur- 
rent SPH simulations (64 3 particles, 
ll/i~ x Mpc box) are probably suffi- 
cient to give accurate predictions of 
the flux decrement distribution at 
these smoothing scales, although a 



factor of eight decrease in the reso- 
lution or box volume would change 
the predictions noticeably. 

Figure ^ compares the flux decre- 
ment distributions predicted by SPH 
simulations of a variety of CDM 
models (with parameters listed in 
Table [lj) to measurements from a 
sample of 28 Keck HIRES spectra. 
The absorption redshifts in the ob- 
servational sample range from z = 
2.5 to z = 3.7, with mean z = 2.96, 
and the total path length in this red- 
shift range is equivalent to 21 full 
Lya to Ly/3 regions. The spectra 
have been obtained by three of us 
(SB, DK, DT) for a variety of pur- 
poses, including studies of the pri- 
mordial deuterium abundance (e.g., 
Buries & Tytler 1998), the statistical 
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Figure 3. The influence of mass resolution and simulation box size on calcu- 
lations of the flux decrement distribution function, in the PM approximation 
at z = 3. Solid lines show results from an LCDM simulation (see Table ^/or 
cosmological parameters) with 64 3 particles in an Uh~ 1 Mpc comoving box. 
In the upper panels, dotted and dashed lines show results from simulations 
with the same box size but 128 3 and 32 3 particles, respectively. In the lower 
panels, dotted and dashed lines show results from simulations with the same 
mass resolution but a 22h _1 Mpc box (128 3 particles) and a 5.5/i — *Mpc box 
(32 :i particles), respectively. 



Table 1. Parameters of the cosmological models shown in Figures ^-|^- By 
definition, Qq — ^cdm + ^ + ^b- 



Model 


| 


If 


h 


OS 





SCDM 


1.0 


0.0 


0.0 


1.00 


0.50 


0.70 


CCDM 


1.0 


0.0 


0.0 


1.00 


0.50 


1.20 


OCDM 


0.4 


0.0 


0.0 


1.00 


0.65 


0.75 


LCDM 


0.4 


0.6 


0.0 


0.93 


0.65 


0.79 


CHDM 


1.0 


0.0 


0.2 


1.00 


0.50 


0.70 



properties of Lyo forest lines (Kirk- 
man & Tytler 1997a), and heavy el- 
ement enrichment in the Lyo for- 
est (Kirkman k. Tytler 1997b). Er- 
ror bars on the data points are es- 
timated by dividing the sample into 
subsets. The statistical uncertainties 
are small because there are many in- 
dependent resolution elements in the 



sample spectra, but the errors are 
not independent. Furthermore, this 
is a comparison between results from 
noiseless simulation data for which 
the true continuum is known and re- 
sults from observed spectra that have 
noise and a locally fitted continuum. 
Substantially more work is required 
to go from this preliminary analysis 
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Figure 4. Comparison of the flux decrement distribution from SPH simula- 
tions of five different cosmological models to measurements from a sample of 
Keck HIRES spectra with mean redshift z = 2.96. Parameters of the models 
are given in Table ^|. The resolution and box size of the simulations corre- 
spond to the solid lines in Figure VI. 



to a definitive comparison between 
models and data. 

Taking the results of Figure ^ 
at face value, there appears to be 
fairly good agreement between the 
observed P(D) and the predictions 
of the SCDM, OCDM, and LCDM 
models. These three models all have 
Gaussian initial conditions and a 
similar mass fluctuation amplitude 
on these scales at this redshift, so 
they predict similar flux decrement 
distributions. There is some discrep- 
ancy with the data for low flux decre- 
ments at 25 km s -1 smoothing, but 
this is the regime where noise will 
have the largest effect on the mea- 
surements. The COBE-normalized, 
Q = 1 model, CCDM, seems to have 
too high a fluctuation amplitude, 
predicting a flux decrement distribu- 
tion that is too broad at 100 km s . 
The cold+hot dark matter model, 
CHDM, seems to have fluctuations 
that are too weak, predicting an ex- 
cessively narrow distribution. 

Another simple continuous field 
statistic is the threshold crossing fre- 
quency jVd, the average number of 
times per unit redshift that the spec- 
trum crosses a specified flux decre- 
ment threshold (Miralda-Escude et 
al. 1996; Croft et al. 1997; Kim ct 
al. 1997). This statistic is analogous 
to the "genus curve" used to char- 
acterize the topology of the galaxy 
density field. Figure ^ uses the PM 
approximation to illustrate the de- 
pendence of jVd on the amplitude 
(top panels) and shape (bottom pan- 



els) of the initial power spectrum 
P(k). We follow Miralda-Escude et 
al.'s (1996) suggestion of plotting JVj 
against F(D), the fraction of the 
spectrum with flux decrement less 
than D, which makes the model pre- 
dictions nearly independent of the 
constant A (equation |lj) that re- 
lates optical depth to overdensity. 
Like the 3-d genus, the threshold 
crossing frequency is sensitive to the 
shape of P(k), with a bluer power 
spectrum leading to choppier Lyo 
forest spectra and hence to higher 
JVd. When the fluctuation amplitude 
is small, the threshold crossing fre- 
quency is independent of the ampli- 
tude of P(k), but at higher fluctu- 
ation amplitudes non-linear merging 
of structures tends to depress JV^, 
again analogous to the behavior of 
the 3-d genus curve. 

Compared to the flux decrement 
distribution, the threshold crossing 
frequency is more sensitive to nu- 
merical simulation parameters, as 
shown by the tests on PM simula- 
tions in Figure ^ These tests sug- 
gest that our current SPH simula- 
tions significantly underestimate 7V<j 
at 25 kms" 1 because of their finite 
mass resolution and slightly overes- 
timate JV(j at 100 km s _1 because of 
their finite box size. 

Figure |?] compares results from 
SPH simulations to measurements 
from the HIRES sample. The CCDM 
model predicts the lowest value of 
because of its high fluctuation am- 
plitude, while the OCDM and (espe- 
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Figure 5. The influence of the amplitude and shape of the primordial power 
spectrum on the threshold crossing frequency. N^, the number of times per 
unit redshift that a spectrum crosses a flux decrement threshold D in the 
downward direction, is plotted against F(D), the fraction of pixels with flux 
decrement less than D. Upper panels show models with an SCDM P(k) shape 
and varying amplitudes; the heavy line is from an SPH simulation with a = 
1.0 and the lighter lines from the PM approximation. Lower panels show 
results from simulations that have the same initial fluctuation amplitude at a 
Gaussian smoothing scale of0.28h~ 1 M.pc (comoving) but have different P(k) 
shapes. 



cially) LCDM models predict higher 
because of the influence of Qo 
and £7a on the ratio H(z)/Ho, which 
determines the relation between co- 
moving h _1 Mpc and the observable 
units of km s" 1 (recall that is 
measured per unit redshift). Because 
of the numerical uncertainties indi- 
cated by Figure |(], we can draw few 
conclusions from the simulation-data 
comparison at present, except to say 
that the model predictions are in the 
right ballpark at 100 km s -1 . The 
different shape of the observed 
curve at 25 km s — 1 may well be an 
effect of noise, which adds spurious 
threshold crossings at low flux decre- 
ments (the leftmost data point is for 
D = 0.04). A useful comparison at 
this scale will require careful investi- 
gation of noise and continuum-fitting 



effects, and perhaps restriction of the 
sample to the data with the highest 
signal-to-noise ratio. 

4 Recovery of the Mass Power 
Spectrum 

Another obvious continuous field 
statistic is the flux power spectrum. 
As one might expect from the anal- 
ogy between equation (^) and "lo- 
cal" models of biased galaxy forma- 
tion (Coles 1993; Fry & Gaztanaga 
1993; Scherrer & Weinberg 1998), 
the shape of the flux power spec- 
trum on large scales is the same as 
the shape of the underlying mass 
power spectrum P(k). On moder- 
ately non-linear scales, it resembles 
the shape of the linear mass P(k) 
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Figure 6. The influence of mass resolution and simulation box size on cal- 
culations of the threshold crossing frequency, in the PM approximation. The 
quantities plotted are the same as those in Figure |^. The simulations (LCDM 
with varying particle numbers and box sizes) are the same as those in Fig- 
ure ^. 
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Figure 7. Comparison of the threshold crossing frequency from SPH simu- 
lations of five different cosmological models to measurements from the Keck 
HIRES sample. 



rather than the non-linear P(k). Fur- 
thermore, if one chooses the parame- 
ter A in order to match the observed 
(D), then the amplitude of the flux 
power spectrum depends only on the 
amplitude of the mass power spec- 
trum. Croft et al. (1998a, hereafter 



CWKH) exploit these facts in a pro- 
cedure that recovers the linear mass 
power spectrum directly from Lya 
forest data. Once the shape of P(k) is 
measured from the flux power spec- 
trum (or, more precisely, from the 
power spectrum of the "Gaussian- 
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ized" flux; see CWKH), the ampli- 
tude is determined using PM simu- 
lations with Gaussian initial condi- 
tions and this P(k) shape, choosing 
the amplitude for which the simula- 
tions reproduce the observed ampli- 
tude of the flux power spectrum. In 
tests on artificial spectra from SPH 
simulations, the CWKH method cor- 
rectly recovers the shape and ampli- 
tude of the true, linear mass power 
spectrum on scales A > l/i _1 Mpc 
(comoving). 

We have recently applied this 
technique to a set of 19 QSO spec- 
tra, obtained between 1987 and 1994, 
originally for the purpose of study- 
ing chemical abundances in damped 
Lya systems (Pettini et al. 1994, 
1997). The resolution of these spec- 
tra ranges from ~ 0.8 — 2.3A (typi- 
cally ~ 1.5A), and the signal-to-noise 
ratio ranges from ~ 10 — 90 (with 
typical S/N > 40). The analysis is 
described in detail by Croft et al. 
(1998b, hereafter CWPHK). 

In gravitational instability mod- 
els with Gaussian initial conditions, 
matching the observed mass function 
of galaxy clusters requires (jgOj^- 5 ss 
0.5 (White, Efstathiou, & Frenk 
1993; Eke, Cole, k, Frenk 1996). Fig- 
ure ra, from Weinberg et al. (1998a), 
compares the CWPHK determina- 
tion of P(k) at 2 = 2.5 to the pre- 
dictions of cluster-normalized models 
with various values of Qo- We adopt 
the P(k) shape parameter T = 0.2 fa- 
vored by studies of large scale galaxy 
clustering. Models with high Ho have 
low erg and predict a P(k) that is too 
low to match the Lya forest results. 
Models with low f2o have high and 
predict a P(k) that is too high. In 
Weinberg et al. (1998a), we formalize 
this argument to obtain constraints 
with 1<t(2(t) uncertainties 



n 



„ .„ +0.12(0.29) /„n 

0.46 „.,„;„,„' open, (3) 



flo = 0.34 



•0.10(0.17) 
+0.13(0.32) 
-0.09(0.16) 



flat, 



for r = 0.2. The difference between 
open and flat models mainly reflects 
the dependence of the linear growth 
factor on The more general best 
fit result is U + 0.2U A = 0.46 + 
1.3(r - 0.2). Higher F yields higher 
Qo because there is less contribution 
to eg from scales beyond that of the 
CWPHK measurement. 

The main assumptions behind this 



Qo constraint are gravitational in- 
stability and Gaussian initial con- 
ditions. With the more restrictive 
assumptions of inflation and CDM, 
we can obtain other constraints on 
cosmological parameters by combin- 
ing the Lya P(k) with COBE-DMR 
normalization. Figure ^| based on 
Phillips et al. (1998), shows power 
spectra of four COBE-normalized 
CDM models with the inflationary 
power spectrum index n chosen so 
that the predicted power spectrum 
matches the CWPHK result (see 
model parameters in Table U). While 
all of these models can simultane- 
ously match COBE-DMR and the 
Lya forest P(k), the high as val- 
ues of the n = 1 models (TCDM2, 
CHDM2) would imply excessively 
massive galaxy clusters at z = 0, 
in agreement with our earlier, more 
general argument. The right hand 
panel of Figure illustrates the con- 
straints from COBE and the Lya for- 
est P(k) alone in the case of LCDM. 
With our fiducial parameter choices, 
this model reproduces the slope and 
amplitude of the CWPHK P(k) al- 
most perfectly. Changing f2o, n, h, or 
fit, in isolation would change the pre- 
dicted slope and amplitude as indi- 
cated. Because the different param- 
eter changes have nearly degenerate 
effects on the Lya P(k), our mea- 
surement imposes a single constraint 
on a combination of the parameters: 

(^) m 2 - 96 (— r 2 >< 

V0.4/ V0.65/ V0.88/ 



n h h 2 
0.02 



1.0 ±0.5. (4) 



As Figure ^| shows, parameter 
choices that reproduce the measured 
amplitude of P(k) also reproduce 
the measured slope. This "coinci- 
dence" represents a significant suc- 
cess of the inflation+CDM scenario: 
the CWPHK measurement confirms 
its generic prediction of a linear 
power spectrum that curves steadily 
towards fc n_4 on small scales. 

We plan to apply this technique 
to a sample of 99 moderate reso- 
lution QSO spectra originally ob- 
tained by three of us (SB, DK, DT) 
for the purpose of identifying candi- 
date systems for deuterium absorp- 
tion studies. The larger sample size 
will greatly reduce the statistical un- 
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log k/(km s- 1 )- 1 log k/(km s" 1 )" 1 



Figure 8. Constraints on f!o obtained by combining the CWPHK measure- 
ment of P(k) at z = 2.5 with the mass function of galaxy clusters at z = 0, 
for Qa = (left) and fl^ = 1 — f2o (right). Filled circles with la error bars 
show the measured P(k). The error bar on the open circle indicates the nor- 
malization uncertainty: at the la level, all points can be shifted coherently up 
or down by this amount. Curves show P(k) at z = 2.5 for cluster-normalized 
models with a power spectrum shape parameter F = 0.2 and various values 
ofQo, as indicated. The cluster normalization constraints, from Eke et al. 
(1996), are listed above each panel. 




-5 -4 -3 -2 -0.8 -0.6 -0.4 -0.2 0.2 

log k/(km s-')-' log A 2 (k p ) 



Figure 9. Constraints on the parameters of CDM models obtained by combin- 
ing the CWPHK P(k) with COBE-DMR normalization. The left panel shows 
linear mass power spectra of four COBE-normalized CDM models (parame- 
ters listed in Table [J), with the inflationary power spectrum index n chosen to 
give the best match to the measured P(k) (filled circles). The right panel illus- 
trates constraints on the parameters of the LCDM model. The asterisk marks 
CWPHK's values of the logarithmic slope 7 and amplitude A 2 = fc 3 P(fc)/2-7r 2 
of P(k) at wavenumber k v = 0.008(km s -1 ) -1 . Concentric circles show the 
la and 2a constraints on these parameters. Our fiducial version of the LCDM 
model predicts the slope and amplitude indicated by the central open circle. 
Filled (open) points show the effect of increasing (decreasing) one parameter 
while keeping the others fixed. Each step corresponds to AQq = 0.1 (circles), 
An = 0.05 (pentagons), Ah = 0.05 (sguares), or Af2b = 0.01 (triangles). 



certainties in P(k), leading to tighter 
constraints on Qq and other param- 
eters. A preliminary analysis of this 
sample yields results within the la 



errors of the CWPHK measurement. 

Figures ^, R ^, and and equa- 
tions (^) and (^) , illustrate the power 
of the Lya forest as a test of cos- 
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Table 2. Parameters of the cosrnological models shown in Figure |^. All of 
these models have Q^h 2 = 0.02. 



Model 


n 


n A 




n 


h 


CT8 




TCDM2 


1.0 


0.0 


0.0 


0.84 


0.50 


0.92 


LCDM2 


0.4 


0.6 


0.0 


0.88 


0.65 


0.74 


OCDM2 


0.55 


0.0 


0.0 


0.88 


0.65 


0.67 


CHDM2 


1.0 


0.0 


0.2 


1.10 


0.50 


0.96 



mological theories. In contrast to the 
current situation for galaxies, for the 
Lya forest we can predict the form 
of the "bias relation" (equation |^) 
based on simple considerations of 
photoionization and adiabatic cool- 
ing. The parameter A that deter- 
mines the strength of this bias is not 
well known a priori, but it can be 
measured for any given model us- 
ing the mean flux decrement, leaving 
all other statistical properties of the 
forest as independent predictions. 
While these statements are based on 
the approximate picture of the Lya 
forest outlined in §2, hydrodynamic 
simulations can provide the neces- 
sary tests of and corrections to this 
picture. The most significant cor- 
rections — peculiar velocities, ther- 
mal broadening, and shock heating 
— are straightforward to compute 
with such simulations. The factors 
that ultimately limit the comparison 
between models and data are likely 
to be observational uncertainties in 
continuum fitting, which make it dif- 
ficult to measure weak fluctuations 
on large scales, and theoretical un- 
certainties in the spatial uniformity 
of the p — T relation. At present it 
appears that these uncertainties are 
small compared to the differences be- 
tween cosrnological models. 

What can we expect to learn from 
the Lya forest over the next few 
years? Most clearly, we should get 
precise constraints on the amplitude 
and logarithmic slope of the linear 
mass power spectrum on scales ~ 
1000 km s" 1 at z ~ 2 - 3.5, from 
application of the CWKH method 
to larger data sets and from inde- 
pendent checks with HIRES data us- 
ing statistics like the flux decrement 
distribution and threshold cross- 
ing frequency. Consistency among 
these statistical properties will test 
the hypothesis of Gaussian initial 
conditions. The threshold crossing 



frequency provides additional con- 
straints on Ho and Qa because of 
their influence on the conversion 
from comoving /t~ 1 Mpc to redshift. 
If continuum fitting can be done with 
sufficient accuracy, it may be possi- 
ble to measure curvature of the linear 
power spectrum or other departures 
from a power law shape. It should 
be possible to detect the signature 
of gravitational growth of fluctua- 
tions over the range z ~ 2 — 3.5. Us- 
ing the rate of growth to constrain 
Qo and Ha probably requires ac- 
curate P(k) measurements down to 
z ~ 1, which will be more difficult to 
achieve because such measurements 
must be based on HST data and be- 
cause the physics of the Lya for- 
est becomes somewhat more compli- 
cated at low redshift (Dave et al. 
1998). Measurements of the flux au- 
tocorrelation or flux power spectrum 
towards QSO pairs can constrain 
through its influence on spacetime 
geometry (CWKH; Hui, Stebbins, & 
Buries 1998; McDonald & Miralda- 
Escude 1998). Finally, all of these 
results can be combined with com- 
plementary constraints from CMB 
anisotropies, the cluster mass func- 
tion, Type la supernovae, and so 
forth to test the inflation+CDM sce- 
nario of structure formation and to 
determine the parameters of the uni- 
verse that we live in. 
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